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Abstract 

It has been suggested recently that the microcanonical entropy of a system may be accurately 
reproduced by including a logarithmic correction to the canonical entropy. In this paper we test 
this claim both analytically and numerically by considering three simple thermodynamic models 
whose energy spectrum may be defined in terms of one quantum number only, as in a non-rotating 
black hole. The first two pertain to collections of noninteracting bosons, with logarithmic and 
power-law spectra. The last is an area ensemble for a black hole with equi-spaced area spectrum. 
In this case, the many-body degeneracy factor can be obtained analytically in a closed form. We 
also show that in this model, the leading term in the entropy is proportional to the horizon area 
A, and the next term is In A with a negative coefficient. 

PACS numbers: 05.30.-d, 03.65. Sq, 04.70.Dy, 04.60.-m, 02.70.Hm 
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INTRODUCTION 



The entropy of a macroscopic black hole is known to be proportional to the area of its 
horizon |]J , in units of the Planck length squared. It has also been shown by several authors 
using a variety of approaches that the leading order correction to this is proportional to the 
logarithm of the area Q] (also see refs. in j^). Recently, a universal form for the (negative) 
coefficient of this logarithmic term has been obtained in ref. j^] by assuming a power-law 
dependence of the area on the mass of the (non-rotating) black hole. For an isolated black 
hole, it is of course appropriate to consider the microcanonical entropy. For a quantum 
system, the microcanonical entropy may be defined uniquely in terms of the degeneracy of 
the state at a given energy, and it has no fluctuation in energy. The many-body degeneracy 
factor, for any nontrivial system, however, is exceedingly difficult to calculate. For this 
reason, it is desirable to approximate the microcanonical entropy by the canonical entropy 
(the leading term), minus a logarithmic term due to fluctuations in the canonical ensemble- 
averaged energy from the equilibrium value. The main objective of this paper is to test 
this formula quantitatively in three solvable models, where the microcanonical entropy can 
also be calculated exactly. The first two of these are systems of noninteracting bosons, and 
not related to black holes. These many-body systems, however, have eigenenergies that 
depend on a single quantum number, similar to a non-rotating quantum black hole. The 
microcanonical entropy of such systems is calculated exactly, and then compared with the 
approximate (canonical) formula. When the logarithmic correction to the canonical entropy 
are included in the canonical expression, its agreement with the microcanonical entropy 
improves markedly. Next, we consider a canonical area ensemble with equi-spaced spectrum 
and distinguishable area components as model for a non-rotating black hole. The area of 
the surface of the event horizon plays a role analogous to the energy |5[. This equi-spaced 
spectrum was first proposed in the early seventies (see ref. pi) and later confirmed by several 
authors using different techniques Q. The model has been studied earlier by Alekseev 
et al. [8] using the grand canonical ensemble. In this paper, we are able to calculate the 
exact microcanonical entropy for this spectrum in a closed form, and obtain the area-law 
for the entropy. We also show that the next term in the entropy is proportional to In A 
with a negative coefficient. This is also verified using the canonical ensemble when the area 
fluctuation is subtracted out. 
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The origin of the black hole entropy, in theories of quantum gravity, is believed to arise 
from the microstates that are generated by a quantum mechanical operator such 
In standard statistical mechanics, when many particles in the mean-field model are trapped 
in a potential well, the microcanonical entropy of the many-body system at a (quantised) 
energy E n is obtained by taking the logarithm of the number of distinct microstates that 
all give the same energy E n . To be more explicit, E n = ^{jvj n ^i e i> where Ni is the 
number of particles with single-particle energy e^. The set {Ni} n denotes a given occupancy 
configuration of single-particle levels that make up a microstate with total energy E n . There 
may be VL(E n ) such distinct microstates for an energy E n , each denoted by a set {Ni}, and 
the microcanonical entropy is then uniqely defined as 

S(E n ) = k B ln(Q(E n )), (1) 

where the Boltzmann constant ks will henceforth be put to unity. Similarly, in models 
of quantum gravity, the (macroscopic) area eigenvalue A n is taken to be coming from the 
elementary components a«, such that A n = Yl{N } n ^i^i, where iVj is the number of elemen- 
tary components with area Oj [8j. Each microstate is specified by a distinct set {Ni}, and 
there may be Q(A) such microstates for a given area A. The microcanonical entropy is then 

S(A n ) = in(n(40). 

ENERGY SPECTRUM DEPENDING ON ONE QUANTUM NUMBER 

Consider a many-body quantum system with eigenenergies E n that are completely spec- 
ified by a single quantum number n, 

E n = f(n) } 71 = 0,1,2,3,.... (2) 

where we assume f(n) to be an arbitrary monotonous function with a differentiable inverse, 
= F(x), such that n = F(E n ). The degeneracy of the states at energy E n is given by 
Q(E n ), a function characterising the quantum spectrum. At this point, we need not assume 
that this many-body system is described in the mean-field picture. The quantum density of 
the system is defined as 

p(E) = J2tt(E n )6(E - E n ) . (3) 
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Using general properties of the delta function, we write 



5(E - E n ) = 6(E - f{n)) = 5(n - F(E))\F'(E)\ , (4) 

where the prime denotes differentiation with respect to the continuous variable E. The 
quantum density of states of the system is then given by 

n=oo 

p(E) = Q(E)\F'(E)\J2S(n-F(E)) . (5) 

n=0 

We know from the Poisson summation formula 

oo oo 

S(E-n)= exp{2t7rk). (6) 

n=— oo k=— oo 

So, if E > 0, we get 

oo oo 

^5{E-n) = 1 + 2^cos(27rA;£) . (7) 

n n=0 k=l 

We thus obtain [9J], from Eq. (JSJ), 

oo 

p(E) = Q(E)\F'(E) | [1 + 2 cos(27rkF(E))] . (8) 

k=l 

We assume the function Q(E) of the continuous variable E to be smooth. The first term 
on the RHS of the above relation is then the smoothly varying part of the density of states, 
while the second part consists of the oscillating components coming from the discreteness 
of the energy levels. For a macroscopic system with large E, the oscillating part may be 
neglected. We then obtain the important relation 

p(E) = n(E)\F'(E)\ , (9) 

where p(E) denotes the averaged smooth density of states. Now we specialize to a system of 
iV noninteracting particles (or in a mean field) constituting the many-body system. Then the 
degeneracy Q(E n ) is just the number of distinct microstates that all have the same energy 
E n , as described in the introduction. The microcanonical entropy is given by Eq. (JTj), which, 
using Eq. (jHJ), may now be expressed as 

S^Hp(E)\F'(E)n, (10) 

where the oscillating part has been dropped. 
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Our next task is to calculate p(E) for a many-body system. This may be obtained by 
considering the canonical partition function of the iV-particle system, and taking its inverse 
Laplace transform using the saddle-point method |3]. The canonical N— particle partition 
function is given by 

POO 

Z((3) = y2n(E n )eM-(3E n ) = / p(E)exp(-{3E)dE , (11) 

where p(E) is defined by Eq. (J3J). Note that we are using the canonical ensemble only 
as a tool to obtain p(E) by Laplace inversion with respect to (3, which is just a variable 
of integration along the imaginary axis. The saddle-point approximation simply gives the 
smooth part p(E). The well-known result is 

p(E) = ^Hl, (12) 

where 

S C (P )=P E + log Z(A>) (13) 

is the canonical entropy evaluated at the stationary point /3 , the prime denotes differentia- 
tion with respect to (3. The energy E is related to the saddle point fi via 



. d\nZ\ 

E = -[^r) • (14) 

Po 



d/3 



Using Eqs. (JTUJ) and (JT2J), we then obtain 



S(E) ~ ScW ~ \ ln(27r^(/3 )) - ln(|F'(^)l)- (15) 

This formula that was originally suggested in Q had missed the last term on the RHS. This 
has also been pointed out earlier in ref. jj]. It turns out that for two of the models that we 
consider in this paper, F'iE) = 1, and the last term in Eq. ()15j) does not contribute. On 
the other hand, in the model with a logarithmic energy spectrum, this term plays a crucial 
role. Note that within the canonical formalism, Sq((3q) = ((E 2 ) — (E) 2 ) is the fluctuation 
squared of the energy jl^ . When this energy fluctuation is subtracted out from the canonical 
entropy, as in Eq. ()15j) . we obtain an estimate for the microcanonical entropy (for quantum 
gravitational fluctuations, see e.g. |llj|). 

The approximation ()15j) for the microcanonical entropy S(E) is very useful, since it is 
prohibitively difficult to calculate it directly from Eq. Generally, in a mean-field model, 



one is given the single-particle quantum spectrum. The direct computation of the many-body 
degeneracy factor Q(E) from this starting point is very time consuming. Instead, it is much 
simpler to obtain the canonical iV-bo dy partition function by well-known recursion relations 
(depending on the quantum statistics) |12| . and then compute the canonical entropy S c (f3). 
Going one step further, one may calculate the canonical energy fluctuation, and use Eq. (JT5J) 
to obtain S. By following this canonical route, no computation of Q(E) is necessary. The 
approximate formula (j!5|) . relevant to black hole physics, has not been yet explicitly tested. 
The main objective of this paper is to test this formula quantitatively in some model systems 
where the microcanonical entropy S can be calculated exactly. 



The logarithmic energy spectrum 

In the first idealized example, we consider iV noninteracting bosons (N — > oo) occupying 
a set of single-particle energy levels (and also the ground state, which is at zero energy) : 

t v = lnp , (16) 

where p runs over all the prime numbers 2,3,5,.... As recently pointed out in ref. 
the many-body microcanonical entropy S of this system is exactly zero. This is because of 
the fundamental theorem of arithmetic, which states that every positive integer n can be 
expressed only in one way as a product of the prime number powers |3| 

n = (17) 

where the p r 's are distinct primes, and n r 's are positive integers, including zero, and need not 
be distinct. It immediately follows from Eq. (JT7j) that the eigenenergies of the many-body 
system are given by 

E n = \nn = ^^n r lnp r , (18) 

r 

and that each eigenstate is non- degenerate. This means that for every macro-state of the 
many-body system, there is exactly one microstate, and Q(E n ) = 1. This implies, by Eq. (flj). 
that the microcanonical entropy S(E n ) = 0. We would now like to check if this is verified 
from Eq. ()15j) . For the many-body spectrum given by Eq. ()18j) . noting that the inverse 
function F(E) = exp(E), we immediately obtain the density of states 

oo 

p{E) = e E [l + 2 cos{27ike E )] . (19) 

k=l 



The second term on the RHS is the intrinsic quantum fluctuation, due to E^s taking only 
discrete values. The smooth part of the density of states is p{E) = e E . To obtain the 
canonical entropy from Eq. (|13|). we need to calculate the canonical partition function Z(f3). 
The exact Z(f3) for N — ► oo in this case is the Riemann zeta function (((3) = Yl^=i n ~^ '> 
which includes the quantum fluctuations. We pick up the smooth part of Z(j3) by evaluating 
it using Eq. (fTTj) with p(E) = e E , obtaining 

Z{p) = {p-\)-\ (3>1. (20) 

From this, we get S c (j3) = — ln(/? — 1) + (3E, so the saddle-point is given by /3 = (1/E + 1). 
Thus the equilibrium canonical entropy is 

S c ((3 ) = E + \nE + l , (21) 

that contains both a linear and a logarithmic term. Evaluation of the fluctuation term is 
elementary, and the microcanonical entropy using Eq. ()15j) is 

S(E) = E + \nE + l-^ ln(2vr£ 2 ) - E = 1 - i 1h(2tt). (22) 

We see that the E dependent terms in the canonical entropy are entirely canceled by the fluc- 
tuation term, the small residual constant is due to the use of the saddle-point method. This 
example is atypical, because the canonical term contains both the linear and the logarithmic 
terms, and still Eq. (|15j) yields (almost) the correct microcanonical estimate. 



The power-law single-particle spectrum 

For our second example, we consider N noninteracting bosons confined in a mean field 
with a single particle spectrum given by e m = m s , where the integer m > 0, and s > 0. 
The energy is measured in dimensionless units. This model is considered here because the 
canonical partition function (for N — > oo) is exactly known Q|: 

oo 1 

ZW = n [l-e^-ffn,)] ' < 23 > 

771=1 L x ' J 

We can therefore calculate the exact microcanonical partition function Q(E) by expanding 
this Z(f3) as in (fTTj) . This is illustrated in the Appendix by taking a quadratic single- 
particle spectrum (s = 2), and showing that the exact combinatorial result for Q(E) of 
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the many-body system can be reproduced by expanding the canonical partition function 
above. It also illustrates the important point that even though the single-particle energies 
are not equi-spaced, the many-body system has equi-spaced eigenenergies E^s. This has the 
consequence that F'(E) = 1, and Q(E) = p(E) when the quantum oscillations are dropped. 
The numerical calculations for a large number of particles, using different power-law single- 
particle spectra, were done in a different context in Ref. ^J, where the details may be 
found. We test the accuracy of Eq. (JT5j) by comparing it with the exact S(E) from Eq. ([TJ 
for s = 1,2. The quantum oscillations have been included in the exact microcanonical 
calculations, but are difficult to see in this scale. In Figs. 1(a) and (b), the dashed curve 
denotes the canonical entropy Sc{E) as given by Eq. (|13|) . and the continuous curve the exact 
microcanonical entropy S (E) for the two power laws. We see from these curves that the two 
differ substantially as a function of the excitation energy E, specially for s — 2. Inclusion 
of the logarithmic correction to the canonical entropy using Eq. (|15j) results, however, in 
almost perfect agreement, as shown by the dot-dashed curves in these figures. 



A MODEL FOR AN AREA ENSEMBLE OF A BLACK HOLE 



The above bosonic model with the power law spectrum e n = n s is not directly applicable 
to the black hole problem since 



3s + l 



p(E) oc E 2 (s+!) exp k(s + 1)E^+ 

where k is independent of E [15I 16], and the leading term of In p(E) does not go linearly 
with E (or A, when taken over to the area ensemble). Following [^j], we consider instead 
distinguishable elementary components, and consider the situation where the elementary area 
components are equi-spaced, aj = j with j taking values 0,1,2,.. etc. For a macroscopic 
black hole with a horizon area A, we assume that the number of independent components 
is 

4 

(24) 



where 7] is a positive constant, and l p the Planck length, which is set to unity. Both iV and A 
are fixed quantities in the microcanonical picture. As shall be shown shortly, the advantage 
of this model is that the expression for the multiplicity Q(A, N) may be explicitly found, and 
therefore the exact microcanonical entropy may be calculated directly. By expanding the 
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microcanonical entropy for large A, we find that the leading term is proportional to the area 
A of the horizon, and the next term goes like In A. This expression is exactly reproduced 
in the canonical ensemble calculation, when the ensemble averaged (A) is identified with 
A, and the fluctuation in (A) is subtracted out. For simplicity, we take the aj-spectrum to 
be nondegenerate. This is identically the one- dimensional harmonic oscillator spectrum in 
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is equivalent to a 



statistical mechanics. The degeneracy g(j) = (j + 1), considered in 
two-dimensional harmonic spectrum and is also examined. For the first case, the one-body 
partition function is 

1 1 



Zi = ^2g(j)exp(-aaj) = ^exp(-aj) 



i — e~" i — x 

3 3 

where a is a variable canonical to the area, and x = exp(-a). The canonical iV-particle 
partition function for distinguishable elementary components is 

Z N = (Zr)" = (1 - x)- N , 



E 



2 3! 



A\ 

A=l 



-X 



J2n {1) {A,N)e- aA , (25) 



.4 



which is analogous to Eq. ([lip. The multiplicity of states of area is therefore 

nV>(A 1 N) = TS **(*! + i) ■ (26) 

This is the microcanonical partition function. The superscript (1) is to distinguish from 
the degenerate case, to be discussed later. It is not difficult to check by combinatorics that 
f^ 1 ) (A, N) for a given A and iV is indeed given by Eq. (|26|). Note that unlike the case of bose 
statistics where the multiplicity Q(E) was found only by expanding the partition function 
(Eq. (JUJ)) or by exact counting, here due to the distinguishability property of the system it 
is given by an explicit formula. Thus an analytical expression for the microcanonical S(E) 
may be found directly from Eq. (|2b1) : 

S{E) = lnfi (1) (A,iV), 

A-l 

= ^ln(iV + f) - In A!, 

N A 

~ Aln(l + — ) +iVm(l + — ) - l/21n(A + A 2 /N) - l/21n(27r), (27) 

J\ iV 



where we have used the Stirling's series and the Euler-Maclaurin |l7] summation formula to 
approximate the sum in step 2 above. 

We now calculate the canonical entropy and show that inclusion of the logarithmic cor- 
rection term (Eq. 1)1 5|l) gives a formula that agrees with Eq. ()27j) for the microcanonical 
entropy. The canonical calculations are performed for a fixed N, and the ensemble averaged 
area is given by 

W = (28) 

at the equilibrium ocq. For large A, we identify (A) = A, and later correct for the fluctuation. 
The canonical entropy is 

S c (a) = aA + In Z N = aA - N\n(l - e~ a ). 

The saddle-point is obtained from the condition that S' c (a ) = 0, and gives 

N 

« = ln(- + l). (29) 

Therefore, 

N 

S 'c^ = (eQ _ 1)2 e " = A i l + A IX)> (30) 

N A 
S c (a ) = Aln(l + -) + AMn(l + -). (31) 

We may now estimate the microcanonical entropy S from Eq. [)15|h where the area fluctua- 
tions have been subtracted from the canonical S c . In the latter case, we get 

S(A,N) ~ 5 c (a )-l/21n^(a )-l/21n(27r), 

N A 

~ Aln(l + — )+A^ln(l + — ) - l/21n(A + A 2 /iV) - l/21n(27r). (32) 
J\ iV 

This is identical to what we obtained in Eq. (|27|) from the asymptotic expansion of the 
microcanonical entropy. Finally, we use the relation (|2*Hl to eliminate iV from the above 
equation. We immediately obtain 

S(A) ~ £A - 1/2 In A - 1/2 1h[2tt(1 + I/77)] , (33) 

where ^ = ln[(l -+- 77) (1 -+- 1/rj) 11 ]. Note that the leading term, proportional to A, comes from 
the canonical entropy S c , and the correction In A arises from the fluctuation in (A). We 
next consider the second case, where the area spectrum is degenerate. The degeneracy is 
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given by g(j) — j + 1. The same case has been dealt with by the authors in Ref |8J, but 
the results were derived in the grand-canonical ensemble. The one-body partition function 
in this case is simply 

z, 1 



(1 -xf ' 

All the calculations including the expression (jHHj) follow as in the nondegenerate case, with 
N replaced by 2N. 

We have shown that for a class of statistical mechanical systems with power-law single- 
particle spectrum, the log-corrected entropy formula (JT5j) accurately reproduces the micro- 
canonical entropy of these systems. Our results are applicable to a large class of black holes, 
with uniformly spaced area spectrum. These include the ones considered in refs. 
Although in the context of loop quantum gravity, the spectrum found in [l8j is not strictly 
uniform, it is effectively equi-spaced for large areas. Moreover, as shown in]Jand [3], exact 
equi-spaced spectrum may emerge in loop gravity as well. It would be interesting to extend 
our analysis to include charged and rotating black holes, which are described by more than 
one quantum number. 
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APPENDIX 

We give an example of the calculation of the exact entropy for a system of A non- 
interacting bosons with a single-particle spectrum given by e m = m 2 , m = 0, 1,2, .... Ini- 
tially, at T = (or E = 0), the particles all reside in the ground state, where m — 0. 
Denote by N ex the number of particles in the excited states. An excitation energy E may 
be shared by N ex out of A particles such that E = Yli N ex = £V ^i- The multiplicity 
Q(E, A) is the number of ways of doing this. For example, take E — 8 and A = 8, then 
there are three distinct configurations. First, N ex = 2 particles can be excited in which case 
each takes 2 2 = 4 quanta and goes to the second level above the ground state. In this case, 
N 2 = 2 and Aj = 0, % ^ 2. Second, N ex = 5 particles can be excited, four of which each takes 
one quantum to the first level (Aq = 4) above the ground state and the other takes 2 2 = 4 
quanta to the second level (A 2 = 1). Finally, all 8 particles can be excited, each takes one 
excitation quantum to the first level (Aq = 8, and Aj = 0, i ^ 1). The energy in all three 
cases are the same: 

E = 8 = Aqei + N 2 e 2 + ... = + 2 x 2 2 + ... + 0, or 

= 4 x l 2 + 1 x 2 2 + ... + 0, or 
= 8 x l 2 + ... + 0. 

Hence, f2(8, 8) = 3. We see that this problem is identical to counting the number of ways 
that an integer E can be partitioned into sum of squares. Note that had we taken 5 < N < 8 
in this example, then there would only be 2 configurations instead of 3, since the last case 
in which each particle takes one quantum is eliminated. In number theory this is known 
as restricted partitioning as opposed to the unrestricted case considered above. Clearly, as 
long as N > E, the number of accessible microstates may be enumerated as if N = oo. 
In Table I below we enumerate the multiplicity for several values E, assuming unrestricted 
partitioning. 

Few remarks are in order. First, as mentioned before, as long as E < N, the enumeration 
of the multiplicity is A-independent. Second, as illustrated above, E is given by a set of 
consecutive integers, E n = n, even though the single-particle energy spectrum is not equi- 
spaced. Each many-body energy level E n has a degeneracy Q(E n , A). This is in fact general 
for any power-law single-particle spectrum and non-interacting particles. Note that the 



12 



E 


N 


uW(E,N ex ,N) 


n(E,N) 


1 


= l 2 


1 


1 


1 


2 


= l 2 + l 2 


2 


1 


1 


3 


= 12 + 12 + X 2 


3 


1 


1 


4 


= 2 2 


1 


1 






= l 2 + 12 + i2 + ]_2 


4 


1 


2 


5 


= l 2 + 2 2 


2 


1 






= 12 + 12 + 12 + X 2 + X 2 


5 


1 


2 


6 


= 12 + 12 + 2 2 


3 


1 






= l 2 + ... + l 2 


6 


1 


2 


7 


= I 2 + l2 + 2 2 


3 


1 






= l 2 + ... + l 2 


7 


1 


2 


8 


= 2 2 + 2 2 


2 








= 12 + X 2 + 12 + X 2 + 2 2 


5 








= l 2 + ... + l 2 


8 
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TABLE I: Calculation of the multiplicity Q(E,N) for iV bosons at an excitation energy E. The 
single-particle energy spectrum is given by e m = 0, 1,4, 9, m 2 . For a given integral E, the parti- 
tioning into the number of N ex excited particles are tabulated in column 1, and the corresponding 
number N ex in column 2. The microstates uj(E, N ex , N) are enumerated in column 3. The last 
column gives the multiplicity Q(E,N) = J2n £X =i ^(E, N ex , N). It is to be noted that for a large 
excitation energy E (not considered in the table), oj(E, N ex , N) may take on values larger than 
unity. 

multiplicity Q(E, N) enumerated in the Table is the same as the expansion coefficient of the 
partition function, i.e. 

TV j 

z ^ = n [i_^]> 

m=l 1 1 

= 1 + x + x 2 + x 3 + 2x 4 + 2x 5 + 2x 6 + 2x 7 + 3x 8 + 4x 9 + 
where x = e _/3 , and the power of x corresponds to the many-body energy E n . 
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